Nearly flat band with Chern number C = 2 on the dice lattice 
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We point out the possibility of a nearly flat band with Chern number C — 2 on the dice lattice in a 
simple nearest-neighbor tightbinding model. This lattice can be naturally formed by three adjacent 
(111) layers of cubic lattice, which may be realized in certain thin films or artificial heterostructures, 
such as the SrTiOa/SrlrOa/SrTiOs trilayer heterostructure grown along the (111) direction. The 
flatness of two bands is protected by the bipartite nature of the lattice. Including the Rashba spin- 
orbit coupling on nearest-neighbor bonds separate the flat bands from the others but maintains their 
flatness. Repulsive interaction will drive spontaneous ferromagnetism on the Kramer pair of the flat 
bands and split them into two nearly flat bands with Chern number C = ±2. We thus propose that 
this may be a route to the quantum anomalous Hall effect and further conjecture that the partial 
filling of the C — 2 band may realize exotic fractional quantum Hall effects. 

PACS numbers: 71.10.Fd,73.43.Cd,73.20.At 
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A few years after the experimental discovery of inte- 
ger quantum hall effect (IQHE) 1 , Haldane wrote down 
a tight-binding model on the honeycomb lattice with 
IQHE 2 , explicitly showing that the essence of IQHE is 
not the external magnetic field. However, it takes more 
than two decades for people to show that the similar 
statement is also true for fractional quantum Hall ef- 
fect (FQHE). Recently several groups have proposed to 
realize FQHE without Landau levels^— . The basic idea 
is to engineer a nearly flat band in two dimcnsions(2D) 
with nonzero Chern number. Electron interaction in this 
partially filled band may realize fractional quantum Hall 
effect, as suggested by exact diagonalization studies^— £. 

In these proposals nearly flat bands are obtained by 
fine-tuning ratios between nearest-neighbor (NN), next- 
nearest-neighbor (NNN), and even further neighbor tight- 
binding parameters. In this paper we point out a route 
to get completely flat bands without this fine-tuning by 
employing a bipartite lattice with unequal number of two 
subsets of sites^. As a concrete example we consider the 
dice lattice as shown in Fig.Q] It is bipartite with unequal 
number of two subsets of sites (the coordination-number- 
3 sites are twice as many as the coordination-number-6 
sites). This system is inversion symmetric with respect 
to the coordination-number-6 sites. We consider a single 
s-orbital with spin-1/2 degrees of freedom(DOF) on ev- 
ery site, and mainly focus on systems close to half-filling, 
i.e., one electron per site. The NN tight-binding model, 
including the Rashba-type spin-orbit coupling(SOC) con- 
sistent with lattice symmetry, will produce two com- 
pletely flat bands separated from the other bands. Be- 
cause the two flat bands are half-filled, ferromagnetism is 
a natural consequence of correlation^—, which gives rise 
to a Zeeman field on the mean-field level. We demon- 
strate the spontaneous ferromagnetism by a variational 
wave function study of Hubbard interactions. 

As a nice feature of the current model system, even 
a small Zeeman field can split this Kramer pair of flat 
bands and produce two separated nearly flat bands with 



Chern number C — ±2. Filling one of them will then pro- 
duce quantized anomalous Hall(QAH) effect with a xy — 
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2^-. This Zeeman field could also be extrinsic, e.g. grow- 
ing the system on a ferromagnetic substrate. Note that 
in a usual ferromagnetic system, a realistic Zeeman split- 
ting would not completely separate the two bands with 
opposite spin polarizations, and a ferromagnetic metal 
results. This is partially why the QAH insulator, which 
needs to be a ferromagnetic insulator, has not been re- 
alized experimentally so far. The main advantage of the 
presently studied system is the existence of the half- filled 
flat bands, which natually support well-separated bands 
by a realistic Zeeman splitting. 

Material realizations: This model hamiltonian may ac- 
tually be relevant to some real systems. Heterostructures 
of transition metal oxide(TMO) perovskites, whose crys- 
tal structures are cubic, are becoming available owing to 
the recent developmem— r— in the fields of oxide super- 
lattices and oxide electronics (for a review, seei£). In par- 
ticular, layered structures of TMO heterstructures can 
now be prepared with atomic precision, thus offering a 
high degree of control over important material properties, 
such as lattice constant, carrier concentration, spin-orbit 
coupling, and correlation strength. 

TMO heterostructures grown along the (111) direction 
have been synthesized experimentally (e.g., Re L 17 ' 18 ). 
Recently it was pointed out that TMO (111) bilayer het- 
erostructures are promising candidates hosting various 
topological phases of matter—. The dice lattice here can 
be formed by three adjacent (111) layers of cubic lattice, 
each of which is a triangular lattice (Fig. [I}. Although 
we considered only simple s-orbital on every site here, 
the result should be valid if the active orbital is a one 
dimensional representation of the D^d group. Some ex- 
amples are the p^-orbital (p x + p y + p z ), and the ai g 
orbital (d yz + d zx + d xy ) of d-electrons under cubic and 
trigonal crystal potential. 

A particularly relevant example is the transition metal 
oxide SrTiOs/SrlrOa/SrTiOs trilayer heterostructure. 
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Note that although the crystal structure of the bulk 
SrlrOa is a monoclinic distortion of the hexagonal 
BaTi03 structure^, thin films of perovskite SrlrOs have 
been synthesized on substrates 2 ^, which are reported to 
be metallic* 2 ^. This indicates that an itinerant electronic 
model could be a good starting point of describing the 
SrTiOs/SrlrOs/SrTiOa trilayer heterostructure. Due to 
the strong spin-orbit coupling on the Ir 4+ ion, together 
with the octahedral crystal field, the active orbital is a 
half- filled effective J e g = 1/2 doublet 21 . The explictly 
form of these doublet in the presence of cubic symmetry 
is \J Z = 1/2) = -±=(+i\xy,}) - \xz,i) +i\yz,\)), and 

\J Z = —1/2) is its time-reversal partner. These half- 
filled orbitals hop around the dice lattice, and contribute 
to states close to the fermi level, indicating the corre- 
lation in the bulk system is intermediate. In a (111) 
heterostructure, cubic symmetry is reduced to trigonal 
symmetry. Nevertheless, to the leading order with re- 
spect to trigonal distortion, the nearest neighbor hop- 
pings between these J = 1/2 orbitals are identical to the 
hoppings of the s-orbitals, which form the model hamil- 
tonian considered in this paper. 

Therefore we think that our proposal is a promising 
route to realize QAH effect. In the same spirit of previous 
works on the FQHE without Landau levels, we conjecture 
that fractional filling of these bands might produce exotic 
fractional quantum hall (FQH) states. The nature of 
these FQH states remains unclear and we leave it as a 
subject of future research. But it is worth pointing out 
that, in a nearly flat band with Chern number C = 2, the 
natural candidate ground states for v — 1/m (m is odd 
integer) filling fractions are non-abelian states described 
by SU(rn)2 Chern-Simons effective theory 2 ^. 

NN model without SOC. The dice lattice and coordi- 
nation system are defined in Fig. [TJ Label the three sub- 
lattices by Vi.,2,3 respectively. Consider a single s-orbital 
with spin-1/2 DOF on every site. As a warmup consider 
NN spin-independent hopping only, 

H ° = - J2 ( tc lc ja + h-c.) - en *> w 

<ij>,a i£\>3 

where a =j, J, labels spin, i, j label sites, rii = J2 a c \a c ia 
is the electron density on site i. Note that sublattice-3 
(V3) has a onsite energy difference e from the other two 
sublattices, as is allowed by symmetry. In this section 
the spin DOF will be omitted. 

In momentum space the hamiltonian reads 

/ i 7 J\ 
ffo(k) = - t lk , (2) 
W k *7k c / 

where 7k = 1 + e ikl + e ik2 , i = *J-i, fc x , 2 = 
k • ei,2 respectively, and the basis is (cik, C2k, C3k)- 
This model has three bands with dispersions E\ — 
- e/2 - x /e 2 /4 + 2t 2 | T i c | 2 , E 2 = 0, and E 3 = -e/2 + 
v / e 2 /4 + 2f 2 |7 k | 2 , as illustrated in Fig. da). 




FIG. 1: (Color online) Top: The dice lattice. Small up- 
ward triangles (bottom layer), downward triangles (top layer), 
and hexagons (middle layer) indicate the three sublattices. 
1,2,3 label the three basis sites in the unit cell at ori- 
gin. Coordination- number-3 sites (1 and 2) and coordination- 
number-6 sites (3) are the two subsets of this bipartite lattice. 
Blue dash arrows labelled by ei and e2 indicate the two trans- 
lations of the dice lattice. Thick green arrows labelled as Dij 
indicate the Rashba SOC directions on those bonds ij, with 
coordination-number-6 site j. Red dotted arrows with the 
labels +x, +y, +z indicate of the projection of cubic lattice 
axis. Captial X, Y, Z are axis for spin space in Rashba SOC. 
Z is the original (111) direction. Bottom: Perspective view of 
three adjacent (111) layers of cubic lattice. The middle layer 
has different color for easy recognition. The top view of this 
tri-layer is the dice lattice. 

The middle band is completely flat as required by 
the bipartiteness. However the top band touches the 
flat bands quadratically at the Brillouin zone corners 
±K = ±(&i = 47r/3,fc2 = 27r/3), similar to double-layer 
graphene^ or certain other models with flat bands 2D . 
The effective two-band hamiltonian at the band touching 
point ±K is 

3t 2 f\Sk\ 2 Sk 2 ± \ 4 . 

where 6~k — k =p K, Sk + = e 17r ^ 3 (Skx + iSky) and 5k- = 

(Sk+y, 

The flat band has Bloch wavefunction (7^,— 7 k ,0) on 
the three sublattices. It has local Wannier functions re- 
siding on the six neighbors of a coordination-number- 
6 site (sublattice-3) with opposite amplitudes between 
sublattice-1 and sublattice-2. 

NN model with Rashba SOC. Rashba SOC induced by 
electric fields can be included as 

#0,SOC =H - t iAc L(6y ' °) a p C 3P + h - C -}> ( 4 ) 

<ij>,a,@ 

where er are spin Pauli matrices, is the unit vector 
along the direction of the cross product x r,-j of elec- 
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FIG. 2: (Color online). Dispersions of NN tightbinding mod- 
els on dice lattice along high symmetry directions. Parame- 
ters used are e — 0.6 1, A = 0.3 1, and B = 0.2 1. The bands 
marked as red in the middle are the (nearly) flat bands. Top- 
left corner is the Brillouin zone with the high symmetry lines 
indicated by dash blue lines, (a). Spin- independent hoppings 
only. (b). Spin-independent hoppings plus Rashba SOC A. 
(c). Spin- independent hoppings plus Rashba SOC A and mag- 
netic field B along Z (111) direction. The Chern number of 
each band is indicated. 
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FIG. 3: (Color online) Local Wannier function of one of the 
fiat bands of NN model on dice lattice (black dotted lines) 
with Rashba SOC A (t = 1 for simplicity). The two compo- 
nent vector (blue) on each coordination-number-3 site indi- 
cates spin-up and spin-down amplitudes on that site. Am- 
plitudes on coordination-number-6 sites vanish, therefore the 
parameter e has no effect. The small hexagon is the "guiding 
center" of this Wannier function. 



trie field Ey and displacement ry for bond ij, A is the 
strength of SOC and is uniform on all NN bonds as re- 
quired by translation and D^d symmetry of the tri-layer. 
The D^d symmetry further restricts the direction of Ey 
to be within the vertical [perpendicular to the (111) lay- 
ers] plane containing bond ij. Therefore are all par- 
allel to the (111) layers. Their directions are illustrated 
in Fig. [TJ In momentum space the hamiltonian reads 
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where 7k± = 1 + e'( fel±27r / 3 ) + e i ( fa±4,r / 3 ), and the ba- 
sis is (cikt,ciki,C2kt,C2k4.,C3kt ; C3ki)- It has three dou- 
bly degenerate bands with dispersions E\ = — e/2 — 

y e 2 /4 + 2t 2 |7k| 2 + A 2 (| 7k .-| 2 + |7^7P), E 2 = 0, E 3 = 
-e/2 + v/e 2 /4 + 2i 2 |7k| 2 + A 2 (| 7k ,-| 2 + |7k,+ | 2 ). If A « 
e, the effective four-band hamiltonian at the original 
band touching point ±K is 

3t 2 f\Sk\ 2 <54^ n i 3V3A* (0 k T \ 

9A 2 

+ — [l 4x4 ± t z ® a z \ + 0(\ 2 Sk) + <3(<5k 4 ), 

(6) 

where Pauli matrix r z acts on the sublattices-1, 2 space. 



The mass term t z ® a z has opposite sign between the two 
band touching points ±K, similar to the Haldane model 2 . 
The dispersions are illustrated in Fig. [^b) . There are 
still two completely flat band dictated by the hamiltonian 
structure. 

The flat bands have Bloch wavefunctions 
(*7krk,iA7k + rk,-*7k r k>-^7k-r kJ 0,0) and its 

Kramer pair, where r k = 7 k 7£ + — 7k 7k-- Therefore 
these flat bands also have local Wannier functions. One 
of the Wannier functions is illustrated in Fig. [3J The 
other Wannier functions can be produced by translation 
and time-reversal. Note that the spin-up component 
of the illustrated Wannier function acquires a phase 
27r/3 under six- fold rotation around its "guiding center", 
similar to the cyclotron orbit in Landau levels except 
that the phase here is twice as large, suggesting Chern 
number C = 2. 

Nearly flat band with Chern number C = 2. The pre- 
vious double-degeneracy of the flat bands is protected by 
time-reversal symmetry. Consider magnetic field effect, 

Hq.soc+b = H ,soc — S^Z • Si (7) 

i 

where Sj = (1/2) J2a,p Aa n a (i c i^ is electron spin, the 
Bohr magneton fig is omitted, and g = 2 (g = —2 for 
the Jeff = 1/2 states of Ir 4+ ion 2 ^) is assumed hereafter. 
The field on sublattice-1, 2 may be different from that on 
sublattice-3. For illustration purpose we draw the band 
structure with uniform field Bi = 0.2 1 along (111) direc- 
tion [Z direction) in Fig. HJc). As expected the Kramer 
pair of flat bands split into two nearly flat bands. Direct 
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computation of Chern numbers shows that they carry 
Chern number C = ±2 [Fig. C^c)]. 

There is a simple physical argument that proves the 
C = ±2 for the two nearly flat bands. Let us turn 
off the Rashba-type spin-orbit coupling A for the mo- 
ment. Note that the sublattice-f ,2 form a honeycomb 
lattice by themselves. We could turn on another artificial 
iXicr ■ (1,1,1) spin-orbit coupling between the second- 
neighbors on these two sublattices only, with the same 
signs as the Kane-Mele model^. In this Ai-only model, 
spin rotation along (111) direction is conserved so that 
we could consider each spin-poloarized subsystem sepa- 
rately. Energy gaps at K and —K, the quadratic band- 
touching points, in Fig. [2fa) are opened by Ai. How- 
ever, it is well-known that the spin-orbit energy gap at 
a quadratic band-touching point transfers Chern number 
one between the two bands. Therefore, two spin-orbit 
gaps at K and —K transfer Chern number two instead 
of one as in the Kane-Mele case. Including a Zeeman field 
B along (111) split the C — ±2 bands, and clearly, the 
resulting bands and their Chern numbers must be very 
similar to those of the A-only model [see Fig. Oc)]. Fix- 
ing a Zeeman field B, it turns out that one can adiabat- 
ically connect the Ai-only model with the A-only model 
by interpolation while keeping all the six bands isolated 
from one another. This adiabatic evolution perserves the 
Chern numbers of each bands. We thus prove the Chern 
numbers in Fig.[2][c). 

Spontaneous ferromagnetism. The flat band is half- 
filled if the entire system is at half-filling. Add onsite 
Hubbard interactions in the hamiltonian, 



-ffint = Hi 



o,soc 



5> 



(8) 



If SOC A = 0, by Lieb's theorem 9 the ground state is fer- 
romagnetic with total spin S = (1/2) [(Ni + N 2 ) - N 3 ] = 
(l/2)iV cc n (iVi j2j 3 is the number of sites on sublattice- 
1, 2, 3 respectively, and equals to the number of unit cells 
N cc \\)- With Rashba A there is no known proof of ferro- 
magnetism. We use a variational (mean field) treatment 
of this problem. 

The ferromagnetic "mean field" hamiltonian is just the 
free fermion hamiltonian with magnetic field i?o,soc+s- 
By inversion symmetry we assume field on sublattices- 
1,2 are the same, B 2 = Bi, but may be different from 
that on sublattice-3, B3. The variational wavefunction 
is the free fermion wavefunction by half-filling this mean 
field hamiltonian. We then evaluate the energy expecta- 
tion value of the Hubbard model H mt and try to min- 
imize it with respect to the variational parameters Bi 



and B3. From preliminary numerical results, the system 
is unstable to spontaneous ferromagnetism for infinites- 
mal repulsive U , consistent with the Stoner criterior>i2, 
however the energy gain is very insensitive to the field 
directions. For e = 0.6 1, A = 0.3 1, U = t, and the field 
directions along (111) (Z direction), the field strength is 
Bt = 0.2440 1 on sublattices- 1, 2 and B 3 = -0.0162 1 on 
sublattice-3. The mean field band structure is very sim- 



0.6 
0.4 
0.2 

S 
-0.2 
-0.4 
-0.6 





J 

















0.6 
0.4 
0.2 

5 

-0.2 
-0.4 
-0.6 



-0.5 



k,/it 



0.5 



FIG. 4: (Color online) Left: The two nearly flat bands (red) 
with Chern numbers C = ±2. Parameters are e = 0.6 1, A = 
0.3 1, Bi = 0.2440 1, and B :i = -0.0162 1. Right: Dispersion 
of a cylinder with 32 unit cell open boundary condition along 
e2 and periodic boundary condition along ei, showing the 
edge states between the nearly flat bands. 



ilar to Fig. [2jc) where a uniform B = 0.2 1 is used. The 
two nearly flat mean field bands are drawn in Fig. |4] The 
three occupied mean field bands have total Chern num- 
ber C = 2 and exhibit anomalous quantum Hall effect. 
The edge state on a cylindrical geometry is also shown 
in Fig. |U 

Conclusion. In this paper we discuss a model 
with spin-orbit coupling on the dice-lattice and the 
correlation physics in it. A transition metal oxide 
SrTi03/SrIr03/SrTi03 trilayer heterostructure grown 
along the (111) direction, where this model may be re- 
alized, is proposed. In this system, two degenerate flat 
bands at half-filling are found. Stoner's instability nat- 
urally leads to ferromagnetism and split the two bands, 
which gives rise two nearly flat bands with Chern number 
±2. This indicate a promising route to realize QAHE. We 
further speculate that further doping into the nearly flat 
Chern bands could lead to FQHE without an external 
magnetic field. We hope these results could encourage 
experimental syntheses and characterization of the ma- 
terial proposed here, as well as future theoretical inves- 
tigations on the nature of the possible FQH states. 

FW thanks the Institute for Advanced Study at Ts- 
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Boston College. 
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